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ABSTRACT 

We report on VLBA astrometry and CXO imaging of PSR J0538+2817 in 
the supernova remnant S147. We measure a parallax distance of 1.471q27 ^P*^ 
along with a high-precision proper motion, giving a transverse velocity Vj_ = 
400^73^ km s~^ A small extended wind nebula is detected around the pulsar; 
the symmetry axis of this structure suggests that the spin axis lies 12° ± 4° 
from the velocity vector (2-D), but the emission is too faint for robust model 
independent statements. The neutron star is hot, consistent with the young 
~ 40 kyr kinematic age. The pulsar progenitor is likely a runaway from a nearby 
cluster, with NGC 1960 (M36) a leading candidate. 

Subject headings: astrometry — pulsars: individual (PSR J0538+2817) — su- 
pernovae: individual (S147) — stars: kinematics 



1. INTRODUCTION 



S147 (G180. 0—1.7) is an optically faint shell-type supernova remnant (SNR) located in 
the direction of th e Galactic anti-center. It is highly filamentary and has a radius of_83' 
dSofue et al.l ll980h. Although age estimates vary widely, from 20 kyr (jSofue et al.l Il980l ) 
to 100 kyr ( IKundu et al.lll980l ). S147 is believed to be one of the oldest SNRs with well- 
defined shell structure. I t has been extens ively stud ied at several wayelengths including 



radio (jKundu et al.l Il980l : ISofue et al.l Il980l ). optical (jLozinskayal Il976l : iKirshner &: Arnold 
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19791 ). UV flPhillips et al.l Il98l[) and X -ray (ISouvageot et al.l ll99Cll ). A recent continuum- 
subtracted Ha image (jPrew et al.ll2005l ) reveals much detailed structure in S147. As shown 
in Fig, m the shell i s mostly circular with a few obvious blowouts to the East and the North. 
Gvaramadzd (120061 ) suggested that a faint blowout to the West and the brighter lobe to 
the East define a bilateral axis passing through the center. Most of the bright filaments 
concentrate on the southern half, where the SNR boundary is also sharper. In contrast, the 
northern half is more diffuse and less well-defined. 



Within the 



jound aries of S147, a 143 ms pulsar PSR J0538+2817 was discovered by 



Anderson et al.l (119961 ) in an undirected pulsar survey using the 305 m Arecibo radio tele- 
scope. The pulsar has a large characteristic age Tc = P/2P = 620 kyr and a dispersion 



measure (^DM)-estimate d distance of 1.2 



J0538+2817 and S147, lAnderson et al. 



ipc. Based on the positions and distances of PSR 
1996h sugges t ed an association between the two. 



From the maximum SNR age of 10^ yr, iRomani fc NeJ (120031) argued that P SR J0538+2817 
had a slow initial spin period. This was supported by iKramer et al.l (120031 ) who obtained a 
rough timing proper motion which indicated that the pulsar is moving from near the SNR 
geometrical center with a kinematic age of only 30 kyr, much shorter than the characteristic 
age. To reconcile the discrep ancy, they also sugg ested a large initial spin period of 139 ms 
for the pulsar. Alternatively, iGvaramadzd (120061 ) proposed that PSR J0538+2817 arose in 
the first supernova from a massive binary, while a second supernova p roduced S147. I n X- 
rays, PSR J0538+2817 has been detected by ROSAT All Sky Survey JSun et al.lll996h and 
HRI imagi ng. A 20 ks Chandra ACIS-S observation discovered extended emission around 



the pulsar (IRomani &: Ngll2003l ). while an XMM- Newton observation (iMcGowan et al.ll2003l ) 
reported the detection of pulsed X-ray emission. 

In summary, S147 with an old nearly circular shell and a central pulsar showing extended 
X-ray emission, provides some unique opportunities to test the details of the pulsar/SNR 
connection. Kinematic measurements, along with X-ray morphology and cooling measure- 
ments can probe the history of the pulsar. 



2. PARALLAX AND PROPER MOTION MEASUREMENT 

PSR J0538-F2817 was observed with the NRAO VLBA over 9 epochs between 2002 and 
2006. The observations were conducted between 1.4 and 1.7 GHz, as a trade-off between 
increasing pulsar flux at lower frequencies and improved resolution as well as reduced iono- 
spheric effects at higher frequencies. At each epoch, 4 spectral windows (IFs) each with 8 
MHz bandwidth are observed simultaneously. To retain visibility phase coherence, repeated 
visits were made to the primary phase calibrator source 133' distant from the pulsar. A VLA 



Fig. 1. — Continuum-subtracted Ha image of S147 (jPrew et al.ll2005l ). The arrow shows the 
pulsar's proper motion direction and points to its current position (see §[2]), with tick marks 
indicating the birth-sites of the p ulsar if born 2 , 40 fc 60 kyr ago. The cross marks the SNR 
geom etrical center sugg ested by iKramer et al.l (120031 ). Solid lines mark a possible bilateral 
axis (lGvaramadzell2006l ) and dashed lines show the direction to several candidate birth-sites 
for the progenitor star. 
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survey of the region around J0538+2817 revealed a compact 8 mJy source only 8' distance 
from the pulsar; the two sources are within the same VLBA primary beam and were observed 
simultaneously. The data were correlated in two passes, first centered on the pulsar field, 
and then a second pass centered on the in-beam calibrator. The signal-to-noise ratio for 
the pulsar was boosted by 'gating' (accepting signal from) the correlator on at the expected 
times of arrival of pulses, using current pulse timing solutions obtained for each epoch from 
ongoing observations at the Jodrell Bank Observatory. Incremental phase calibration de- 
rived from observations of the in-beam source were applied to the pulsar for considerably 
improved calibration. 

Data analysis was performed using the AIPS package with a customized pipeline (Brisken 
2007, in preparation), which included amplitude calibration based on the system temper- 
atures and antenna gains, followed by visibility phase calibration based on the primary 
calibrator and the in-beam source. Large scale ionospheric phase effects were estimated 
and corrected with the AIPS task TECOR, using global models of the ionospheric electron 
content based on distributed GPS measurements. Once the calibration was completed, the 
pulsar position was measured by fitting a Gaussian ellipse to calibrated images made sepa- 
rately for each IF. The pulsar parallax and proper motion were then obtained from a linear 
least-squares fit to the epoch positions. 

Table [1] shows the best-fit astrometric results from 8 good epochs, for which the pulsar is 
far away from the sun, with IFs equally weighted. The derived pulsar distance and transverse 
velocity are d = lAI^^ kpc and V±_ = 400!l^73'' respectively. The distance is in good 

agreement with the D M-estimate value c^dm = 1-2 ± 0.2 kpc, using the NE2001 model 
( ICordes fc Lazid 120021). The proper mo tion is also consistent with the timing observation 
results reported by iKramer et al.l ( 20031 ). but our errors are much smaller, especially in the 
ecliptic latitude (/i/3). 

To convert the proper motion to its local standard of rest ( LSR), we correc t for dif- 
ferential Galactic rotation using 

fio = 220 km s"\ Rq = 8.5 kpc flFich et al.lll998h a nd the 
solar constants 10 ± 0.36, 5.25 ± 0.62, 7.17 ± 0.38 km s"^ JPehnen fc BinnevI Il998[ ). The 
corrected proper motion is /i* = 58.51 ± 0.18 mas yr^^ at position angle (PA) 336?8 ± 0?13, 
which converts to a transverse velocity of V± = 4071^4^ ^™ Backward extrapolation 
of the pul sar's motion indicate s that it passed ~ 8' from the geometrical center of S147 
defined bv iKramer eTaP J^OoJ at a = 05M0'^01'^ ± 2^ 6 = 27°48'09" ± 20" (J2000). Since 
the pulsar is young and located near the Galactic plane, this trajectory is not significantly 
altered by acceler ation in the Galactic potential. For any reasonable Galactic model (e.g. in 
Sun fc Hanll2004l . and references therein), the displacement from a linear trajectory over the 
pulsar characteristic age is much smaller than 1'. However, the SNR and pulsar are almost 
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certainly associated (see below) and thus we conclude that the nominal center of the shell 
does not represent the explosion site and that S147 underwent asymmetric expansion, likely 
due to inhomogeneities in the surrounding ISM. 



3. CXO OBSERVATIONS 

Chandra observations of PSR J0538+2817 were carried out on 2006 Jan 18 & 20 (Ob- 
servation IDs 6242 & 5538) with the ACIS-I array operating in very faint timed exposure 
(VF TE) imaging mode. The pulsar was positioned near the aim point on the 13 chip, which 
was the only chip active during the observation. Subarray mode with 160 rows was used to 
further reduce the CCD pile-up. The resulting frame time of 0.5 s reduced the pile-up of the 
pulsar to 2%, ensuring the distortion of the PSF and the high energy spectrum is minimal. 
The total live time was 93.2 ks and examination of the background light curve showed no 
strong flares during the observation. Hence, all data is included in the analysis. 

For comparison, we have also reprocessed the archival 20 ks ACIS-S exposure (ObsID 
2796) observed on 2002 Feb 7. After filtering out the periods suffering from background 
flares, 18.5 ks of clean exposure remains. All data analysis was performed using CIAO 3.3 
and CALDB 3.2.1 to ensure the latest time-dependent gain calibration and charge transfer 
inefficiency (CTI) correction are applied. To further improve the s patial resolution, we 
removed the ACIS pixel randomization and applied the algorithm by 
correct the position of split pixel events. 



Mori et all (|2001[ ) to 



SPATIAL ANALYSIS 



The ACIS-I 0.5-8 keV image is shown in Fig[21 The point source produces 0.12 cts s^^. 
With the short frame time, a trail of counts can be seen along the readout direction (at 
PA 5°), with .5 cts pixel"^. The extended emission ~ 5" NE of the pulsar reported by 
Romani fc Ngj (120031 ) is clearly detected in the new observation. However, the deeper expo- 
sure shows that the structure is not an obvious torus. It also appears somewhat fainter in 
the new data, with (4.3 ± 0.7) x 10"'' cts s~^ in 0.5-8 keV after background subtraction, as 
compared to (9.7 ± 2.2) x 10~^ cts s~ \ in the archiva l ACIS -S data. Diffuse emission is also 



seen in the SE direction, as noted by iRomani fc Ng (120031 ). at ~ 3" from the point source. 
Given that the direction is behind the pulsar's proper motion, this could be trailed emission 
from relativistic electrons in the mot ion- confined pulsar wind, as observed in other PWN 
systems. 
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Fig. 2. — ACIS-I 0.5-8 keV image. The arrow indicates the pulsar's proper motion direction 
in its LSR. The dotted hnes separate the four quadrants used to measure the azimuthal 
distribution of the extended emission. 



Fig. [3] compares the observed surface brightness of the source for the ACIS-I data and 
the model PSF-|-background in 3 energy bands. The PSF is simulated using the Chandra 
Ray- Tracer (ChaRT) and MARX software using the best-fit pulsar spectrum from a 1" 
radius aperture, in order to minimize any nebular contamination. Although the readout 
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Radius (arcsec) 



Fig. 3. — Observed source (points) and model point source (line) radial count distributions 
in three energy bands, with the corresponding uncertainties (shaded bands). 

trail is simulated in the model PSF, we excluded two rectangular regions of I'.'S wide along 
the readout trail and beyond 2'.'5 from the point source in our analysis, in order to improve 
the statistics. The graph suggests that the data and PSF model are well matched for the 
low energy band, while excess counts appear at several radii for the higher energies. 
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To investigate the azimuthal distribution of these counts, we show in Fig. H] the PSFs 
and observed counts per unit area for the four quadrants (in three energy bands) with Quad 
1 along the proper motion direction (see Fig. [2]). The shaded region shows the Poisson 
uncertainty in the surface brightness measurements. Here we see significant departures from 
the PSF in all three energy bands. The excesses from 4" — 6" a nd 9^^ — 11^^ in Quad 2 in 
the low energy band represents the candidate 'torus' suggested in lRomani fc Ngl (120031 ): no 
corresponding excess is seen in Quad 4. A persistent excess is also seen from 4" — 8" in Quad 
3 of the medium energy band, representing the diffuse emission trailing behind the pulsar. 



0.5-1 .5 keV 1 .5-4 keV 4-8 keV 




r (arcsec) r (arcsec) r (arcsec) r (arcsec) r (arcsec) r (arcsec) 



Fig. 4. — Surface brightness of the observed source (points) and model point source (lines) 
for the four quadrants in three energy bands. The y-scale on the left applies to the 0.5-1.5 
keV panel only, while the right scale applies to the other two panels. 



Most intriguing, however, is the very significant excess seen in the 4-8 keV band from 
1" — 4" (Fig. [3]). In Fig. H] we see that this lies in quadrants 2 & 4, indicating significant 
structure with a hard spectrum very close to the pulsar. The 4-8 keV ACIS-I image (Fig. [5]) 
in fact shows a clear symmetric, almost linear structure extending ~ 2'.'5 from the pulsar. 
Realizations of the 4-8 keV model PSF are nearly circular and show no obvious spikes. 
With only 32 counts, we cannot resolve any details or assign a clear origin for this emission; 
however, it should be noted that the X-ray images of many other young neutron stars show 
polar jets or equatorial tori. // this is the case for PSR J0538-I-2817, the symmetry axis 
could indicate the pulsar spin axis. Mor e importantly, its alignment with the proper motion 
could constrain the pulsar kick physics (jRomanill2005l : Ng & Romani 2007, in preparation). 
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In order to measure the PA of the symmetry axis quantitatively, the photon positions were 
fitted to a straight hne passing through the point source. Linear least squares fitting was 
employed and we obtained the best-fit PA at 79° ± 4°. 




Fig. 5. — Left: ACIS-I 4-8 keV image with the dashed line showing the best-fit PA of the 
counts. The arrow indicated the pulsar's proper motion in its LSR. The circle is 3" in radius. 
Right: Model PSF in the same energy band. 



Although the 4-8 keV excess is statistically significant, interpretation of such faint com- 
pact structure is of course very uncertain. Some infor mation on the pulsar orientation may 



be gleaned from the radio polarization measurements of iKramer et al.l (l2003l ). which suggest 
a magnetic inclination a = 95° and impact angle a = 2°, implying a spin axis inclination of 
( = a + a ^ 97°, i.e. very nearly in the plane of the sky. This means that both the polar jet 
interpretation (with the near-orthogonal view implying similar jet /counter-jet fluxes) and 
the equatorial torus interpretation (with a thin nearly edge-on torus) remain viable. If the 
extended emission is interpreted as polar jets, then the pulsar spin axis is 78° off the proper 
motion direction. This is much larger than any of the other PSR/PWN systems. If the axis 
is the equatorial (torus) plane, then the inferred spin vector is at PA 169° (—11°), i.e. 12° off 
the velocity vector. 

A statistical argument in fact supports the latter interpretation. If the spin axis and 
proper motion axis are orthogonal in true space, then most observed orientations make the 
projected angle on the plane of the sky smaller. If they are aligned, the projected angle tends 
to remain nearly aligned. For completely random orientations, 90% of vectors separated by 
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0°-10° retain a 2-D angle < 10°; for vectors separated by 80° -90° only 29% retain 6*20 > 80°. 
Thus the observed angle is, a priori, ~ 3 x more likely to arise from an aligned system. One 
additional piece of evidence could be extracted fro m absolute position angle measurements 
of the radio polarization (c.f. I Johnston et al.ll2005l ). Although mode ambiguity allows a 90° 
jump in the inferred projected field, these authors note that most pulsars appear to emit 
in the orthogonal mode, which may allow one to discriminate between the jet and torus 
interpretation. 



5. SPECTRAL ANALYSIS 

We now turn to the spectral analysis of the point source. For the best possible con- 
straints on the spectrum, we have reprocessed both the 18.5 ks cleaned ACIS-S data set 
and our new 93 ks ACIS-I data with the latest time-dependent gain calibration and charge 
transfer inefficiency (CTI) correction. The source spectrum was extracted from a 2" radius 
aperture with the script psextract, and the response matrix files (RMFs) were replaced by 
the ones built using the tool mkacisrmf , which accounts for the CTI. To model the aperture 
corrections, 10 PSFs with monochromatic energies from 0.5 to 9.5 keV were simulated using 
ChaRT. The fractional energy encircled by the aperture as a function of energy is obtained, 
and then used to correct the ancilla ry response fi les (ARFs). As the ACIS-S data suffer 20% 



pile-up, the CCD pile-up model by iDavid (1200 ll ) is used in all the spectral fits 



Results from the combined fits of the ACIS-S and ACIS-I datasets are listed in Table [2l 
All fits are to the 0.3-8 keV range and the spectral parameter errors reported are projected 
multidimensional 1 — cr values. For the uncertainties in fiux (and hence stellar radius) , as is 
often the case with low-statistic CCD-quality data, the projected errors are very large due 
to uncertainties in spectral parameters. Therefore we followed other authors in reporting 
the single parameter (i.e. 1-D) 1 — a error for the flux. 

The source spectrum is adequately fitted by an absorbed blackbody. The best-fit A^h = 
2.47 X 10^^ cm~^ is lower than the previous CXO results, but consistent with the XMM- 
Newton measurements. Comparison with the DM value of 39.7 pc cm ~^ gives n-a/ne = 24, 
which is relatively large. At the pulsar distance of 1.47 kpc, the best-fit spectral parameters 
give an effective blackbody radius of = 2.19 km. This is too small to be reconciled with 
the whole stellar surface, but the flux could be hot T ~ 2 x 10® K emission from a small 
fraction of the stellar surface (~ 2.7% for an i?oo = 13.1 km star), possibly due to some 
heating mechanism such as bombardment of the polar cap regions by relativistic particles 
from the magnetosphere. The thermal radiation from the neutron star surface could also 
be described by atmospheric models. Light-element neutron star atmosphere models, such 
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Table 1. PARAMETERS FOR PSR J0538+2817 



Parameter 



Value 



Epoch (MJD) 
R.A. (J2000) 
Decl. (J2000) 
Ha cos 5 (mas yr" 
fj.s (mas yr"^) 
TT (mas) 
d (kpc) 
V± (km s-i) 
I, b 



III, lib (mas yr 



-1^ 



Up cm 



53258 

05^38°'25!05237 ± O'^OOOOl 
+28°17'09'.'3030±0'.'0001 
-23.53 ±0.16 
52.59 ±0.13 
0.68 ±0.15 

n 47+0.42 

4001^^" 

179?7186, -1?6859 
-57.03, 8.18 
0.027 ±0.006 



Table 2. SPECTRAL FITS TO PSR J0538±2817 



BLACKBODY/ ATMOSPHERE 



MODEL 


(1021 cm-2) 


(10« K) 


Rao 

(km) 


Abs. Flux 

/0.5-8 
(lO^i-^ ergs 
cm^^g-i) 


Unabs. Flux 

/0.5-8 
(10^1'^ ergs 
cm^^g-i) 


BB 

Atm 
BB+PL 
Atm+PL 


2 47+0.15 
2 Q4+0.07 

2.50 ±0.15 
2.95 ±0.07 


9 -,1+0.03 
^•^^-0.04 

1.05 ± 0.05 

2 10+0.04 

l.UO_Q Qg 


2.19 ±0.01 
11.16 ±0.02 
2.23 ±0.01 
10.99 ± 0.02 


7.24 ±0.07 
7.21 ± 0.08 
7.15 ±0.07 
7.14 ±0.08 


17.4 ±0.2 

21.4 ±0.2 

17.5 ±0.2 
21.4 ±0.2 



POWER LAW (F = 1.5) 

Abs. Flux Unabs. Flux 

/0.5-8 /0.5-8 
(lO^i-^ ergs (10"" 



ergs 



ergs 



x'/doi 



0.20 ± 0.06 
0.20 ± 0.05 



0.24 ± 0.07 
0.24 ± 0.07 



109.7/139 
120.7/139 
99.5/138 
110.2/138 
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as those dominated by hydrogen, have large Wien excesses. This gives a lower effective 
temperature and hence larger stellar radius in the fit. We use here a pure H model with 
B = 10^^ G ( iPavlov et al.lll992l : IZavlin et al.lll996f). as the infer red surface magnetic field 
strength of PSR J0538+2817 is 7 x 10^^ G f lAnderson et al.lll996[ ). During the fit, the mass 
of the neutron star is held fixed at M = IAMq and the normalization constant is fixed using 
the pulsar distance. The best-fit surface temperature and radius are = 1.05 x 10^ K and 
Roo = 11.2 km. As expected, this model suggests a lower effective temperature covering a 
large fraction of the neutron star surface for a canonical radius. The fit is statistically slightly 
worse than that of the blackbody model, but both produce quite acceptable values. 

The extended emission observed at high energies suggests some flux in the 2" aperture 
is contributed by non-thermal emission. Therefore we tried adding a power law component 
with fixed F = 1.5 to the models. Although the results are not improved substantially, 
the non-thermal flux is detected at > 3 — cr level as shown in the table; the effect on the 
parameters fitted for the thermal component is very small. The best-fit atmosphere+power 
law model is shown in Fig. [6l 



6. DISCUSSION 



6.1. PSR J0538+2817 / S147 Association 



The association between PSR J0538+2817 and S147 suggested by I Anderson et al.l (119961 ) 
was based on their positional coincidence and apparent consistency of the distances and 
ages. In particular, the authors argued that since the pulsar location is near the SNR 
center, it is unlikely to be a chance association. With our accurate proper motion mea- 
surement we can improve this argument. In order to have a quantitative estimate of the 
chance alignment probability, w e did simple Monte Carlo simulations using the model by 
Faucher-Giguere fc Kaspil (120061 ). Following these authors, we assur ne the pulsars are born in 
the Galactic plane with the galactocentric radial distribution from lYusifov &: Kiiciikl (120041 ) 
and exponential distribution in the scale height, and with birth velocities distributed as 
a two-component Gaussian model. Acceleration due to Galactic potential was ignored for 
simplicity. Our results show that in 2 x 10^ yr, the most extreme age estimate for S147, 
only 1 pulsar in 2 x 10^ would have a chance passage within 8' of t he S147 cent er. We 



also applied the Maxwellian pulsar velocity distribution suggested by iHobbs et al.l (120051 ) . 
obtaining nearly identical result . With a Galactic neutron star birthrate of 2.8 per century 
(IFaucher-Giguere fc Kaspil l2006l ) and a radio beaming factor of ~ 1/5, the probability of 
finding a random, unassociated radio pulsar younger than 1 Myr which has passed within 8' 
of the SNR center is ~ 3 x 10~^. This estimate is a conservative upper limit to the probabil- 
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12 3 4 

E(keV) 

Fig. 6. — Point-source spectrum with a pileup-corrected magnetic H model atmosphere + 
power law spectrum and residuals. 

ity: we believe the true age of S147 is considerably younger, and the high X-ray temperature 
of PSR J0538-I-2817 also implies a younger < 10^ yr age. Thus a more realistic chance prob- 
ability is > lOx smaller. To conclude, PSR J0538-I-2817 is almost certainly associated with 
S147 and this implies a SNR distance of ~ 1.5 kpc. This value is substantially larger than 
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some previous estimates (e.g. lKundu et al.lllQSOl). thus it calls into question some papers that 



assume a much closer distance to S147 (e.g. iPhillips et al.lll98ll : ISallmen fc Welshll2004j ) 



6.2. S147 as a Cavity Explosion 

In standard SNR evolution, th e shell radius in the Sedov- Taylor phase is given by 
-RsNR = 0.31(i?5i/r2o)^/^t^/^ pc (e.g. Ivan der Swaluwl I2OOII ). where the explosion energy is 
Eq = lO^^E^i ergs, age in t years and externa l medium of density uq cm~^. The observed 
angular size of S147 9 = 83' fISofue et al.lll98Cl) co rresponds to a physical radius of 35 pc at 



1.47 kpc. For an age of 30 kyr (IKramer et al.ll2003l ). this requires a very energetic explosion of 
E^i = 20 no. This suggests S147 probably occurr""' " — j u„uui„ i;u„u. 



Kramer et al. 


2003; 


Gvaramadze 


2006) 



Hence the SNR had a long free-expansion phase, only passing to the Sedov- Taylor phase when 
it reached the cavity boundaries at relatively large radius. This scena rio receives further 



supp ort from the observed low expansion velocity of S147 at 80 km s ^ (IKirshner fc Arnold 



19791 ). Note that the progenitor's proper motion c an make the wind bubble asymmetric. 



with the cavity extending further behind the star (iGvaramadzd l2006l ) . Density gradients 
in the external medium can also enhance this asymmetry. Indeed gas and dust surveys 
suggest that the medium to the South of S147 is denser and the shell is flat tened with 



the br ightest filaments on this side. Supporting this, the optical observations by iLozinskaya 



( 1l976l ) provides some hint that the expansion rate is faster in the northern half of S147. Thus 



we generally expect the 'geometrical center' to lie somewhat North of the true explosion site. 



6.3. Birth-site of the Progenitor 

With the inferred distance to S147, it is possible to search for the birth-site of its parent 
star. O and B stars are the direct progenitors of neutron stars and the minimum mass for 
supernova explosion is ~ 8Mq. These massive stars are generally formed in OB associations 
and young open clusters. With their short < 50 Myr lifetimes they do not travel far from 
the birth-sites. Typical peculiar velocities are a few km s~^; for 10 km s~^, we expect the 
progenitor to travel < 500 pc. We compiled a li st of open clusters and OB a ssociations 
from the catalogs including iRuprecht et al.l (119831 ) , iMel'nik &: Efremovl (119951 ) , iDias et al. 



(2002) & Kharchenko et al. (j2005 ) and found only four candidates younger than 50 Myr 



with a nominal distance to S147 of < 500 pc (Table [3]). Directions to potential birth-sites 
are shown in Fig. [H 
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With a 3-D spatial separation of only 220 pc from S147, N GC 1960 (M36) is the cl osest 
candidate. This is a relatively massive ~ 40 Myr-old cluster (IKharchenko et al.ll2005l ) and 
we consider it the prime candidate birth-site. If correct, the progenitor traveled from this 
cluster at PA ~ 175° with a 3-D space velocity of ~ 5(M/1OM0)^-^ km s~^. We see no 
evidence in gas maps or the SNR shell that the progenitor wind has disturbed the denser 
cloud to the South. This and the blowouts of the SNR to the North argue against Gem 1 
or NGC 2175 as the parent cluster. Stock 8 to the North remains viable, but is distant at 
> 400 pc and less massive. Finally, it is intriguing to note that the 5*^ nearest young cluster, 
the very poorly studied NGC 2331, lies precisely in the direction of the largest (Eastern) 
extension of the SNR shell; the Ha shell here extends 40% (~ 35') further than the main 
shock front, with wispy emission present up to 10' further in this direction. It is tempting 
to associate this blowout with a stellar wind trail extending along the path to NGC 2331, 
but this cluster appears to be relatively low mass and may be too old, with no remaining B 
stars. 



6.4. Explosion Site 

From the precise pulsar proper motion, we know that the supernova explosion must 
have occurred along the line in Fig. [T] with the arrowhead at the present pulsar position. 
The explosion site is determined by the true age of S147/PSR J0538+2817; tick marks 
on the line indicate 20, 40 and 60kyr ages. The simplest interpretation is to infer birth 
at the closest approach to the geometrical center 36kyr ago and assign an uncertainty of 
~ 8'/57.6 mas yr~^ = 8kyr. However, if we can define a second axis for the progenitor 
motion, we can ob tain a more precise age. If we adopt the symmetry axis suggested by 



Gvaramadzd (120061 . solid lines in Fig. 1) then the birth-site is to the NW of the geometrical 



center and the intersection with the proper motion suggests a SNR age near 20 kyr. Given 
Table 3. Open Clusters younger than 50Myr within 500 pc of S147. 



Name 


d (kpc) 


en 


r (pc) 


N* 


Age(Myr) 


NGC 1960 (M36) 


1.32 


6.4 


217 


60 


42 


OB Gem 1 


1.34 


9.9 


276 




< 5 


NGC 2175 


1.63 


10.0 


311 


60 


9 


Stock 8 


1.82 


7.1 


405 


40 


41 


NGC 2331 


1.33 


19.3 


488 


30 


? 
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the rather irregular nature of the Northern half of the remnant and the argument that the 
explosion should be south of the geometrical center, we do not find this axis convincing. The 
axis to NGC 2331 intersects at a more plausible explosion age of 30 kyr. 

If however we adopt NGC 1960 as the birth-site, then a explosion somewhat south of (in 
front of) the geometrical center becomes natural. Without a blowout identifying an entry 
site, we cannot set a precise axis, but the path should pass close to the geometric center, 
implying an intersection with the proper motion vector at < 60 kyr. Our conclusion is that 
the best estimate of the SNR age is 40 k yr with a maximum p lausible range of 20 — 60 kyr. 



This is slightly older than the estimate of iKramer et al.l (120031 ). However, the inferred initial 



spin period of PSR J0538-I-2817 is not significantly changed. Assuming magnetic spin-down 
with constant braking index n = 3, the initial spin period is given by 



n — 1 T 

Pn = P 1 ^ 



1 



-1 



where r and Tc are the kinematic and characteristic ages of the pulsar respectively. We 
obtained Pq = 138 ± 2.3 ms; the kinematic age of the pulsar is indeed much smaller than its 
spin-down age. 



6.5. Spin-velocity Alignment 

Only a few pulsars have estimated initial spin periods; the value for PSR J0538+2817 
is the longest among these. Of course, the high space velocity of the pulsar argues for a 
strong birth kick. Note th at with the large pa rallax distance to the pulsar, the binary break- 



up scenario described by (iGvaramadzd l2006l ) is now even more improbable. If one further 



accepts a progenitor origin in the Galactic plane near NGC 1960, then the pulsar's present 
rapid return to the plane further supports a birth kick uncorrelated with its parent's motion. 

These considerations make a comparison of this pulsar's kick and spin direction partic- 
ularly appealing. Unfortunately the morphology of the extended structure near the pulsar 
is not clear enough to define a definitive spin axis and thus weakens this system's ability 
to cleanly test the spin-kick models. However, we can turn the question around: if the 
mechanism that seems to cause spin and kick alignment in other young pulsars acts on PSR 
J0538-f2817, do we expect its spin to be more nearly aligned or orthogonal? To retain a 
slow spin with a large kick velocity, the net kick vector must be nearly radial, applying little 
torque to the star. 



We have performed a series of simulation of neutron star birth kicks in a range of 
models where a single thrust is applied to the surface at fixed angle as the proto-neutron 
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star cools, with amplitude proportional to the driving neutrino luminosity (Ng & Romani 
2007, in preparation). Comparing with the set of all neutron stars with initial spin and/or 
speed measurements, we have found (for several models of proto-neutron star evolution and 
neutrino cooling) the best-fit parameter distributions for the neutron star pre-kick spin, the 
kick amplitude, normal direction and duration that reproduce the observed pulsar spin and 
speed distributions. For these parameters (fixed by a set of ~ 50 other pulsars) we can ask 
whether a pulsar with slow initial spin like PSR J0538-I-2817 is more likely to have its birth 
velocity and spin vectors aligned or orthogonal. The simulations find that for P ~ 140 ms, 
the aligned case is produced 30 — 90x more frequently than the orthogonal cases, even though 
faster spin pulsars do not always show good alignment. Ng & Romani (2006, in preparation) 
discuss the significance of this result for improving kick constraints. 



6.6. Pulsar Thermal Emission 



From the X-ray spectral results, the neutron star atmosphere fit gives an effective surface 
temperature of > 10^ K. T his matches well to the standard cooling curve for the pulsar age 
of S. 40 kyr (c.f. F i g 6 in iMcGowan et al.ll2003l ). Comparing with the cooling models of 
Yakovlev fc Pethickl (120041 ) we see that our best-fit 1.05 x 10^ K surface agrees well with a 
typical cooling model at age 25 kyr. Even so-called slow cooling neutron stars (low mass 
stars, with crustal neutron pairing and/or accreted low Z envelopes) drop very rapidly below 
Too = 10^ K after 10^ yr. Thus if we interpret the thermal emission as full surface emission, 
it seems impossible for PSR J0538+2817 to be as old as its characteristic age. In contrast, its 
thermal surface emission is quite consistent with its young ~ 30 kyr kinematic age, requiring 
no direct Urea process or any other exotic cooling mechanisms. 

XMM observation s find a low 18% soft X-ray pulsation with a very broad profile 



( IMcGowan et al.ll2003l ) . These authors interpret this as a hot polar cap from a nearly aligned 
rotator (polar cap axis and rotation axis both close to the line of sight). This seems at odds 
with the radio polarization and PWN data, so a more natural interpretation might be emis- 
sion from a gradual temperature variation across a ligh t element surface, perhaps caused by 
magnetic dipole variation in the thermal conductivity (jCreenstein fc Hartkdll983l ). 



If interpreted as a re-heated cap emission the thermal fiux would be a surprisingly large 
1% of the full spin-down power. On the other hand, the non-thermal emission from the 
PSR/PWN system as a whole is close to that expected. In an aperture of radius 15", 
the observed count rate in the 2 — 10 keV band is 5 x 10^^ cts s~^. After subtract- 
ing the background and thermal emission from the pulsar, the non-thermal contribution 
gives 2.8 X 10"'^ ctss~^. For a power law of P = 1.5, this converts to the luminosity of 



- 18 - 



10 erg s . iPossenti et al.l (120021 ) found an empirical relation between the X-ray flux 
in the 2-10 keV band and spin-down luminosity: logLx,(2-io) = 1-34 log Lsd — 15.34. For 
the case of PSR J0538 -I-2817, the spin-down luminosity derived from the radio parameters 
([Anderson et al.l Il996l ) is = 5 x 10^^ ergs s~^, which predicts an X-ray luminosity of 



^X,(2-10) 



10 erg s , very close to the observed value. 



6.7. Conclusions 



We have reported VLBA astrometric measurements and Chandra ACIS-I observation 
of PSR J0538-I-2817. The VLBA astrometry gives the pulsar distance of 1.471q27 ^P*^ with 
a precise model independent transverse velocity V± = 400^^3^ km s~^. These observations 
strengthen the association with S147 and suggest NGC 1960 as plausible birth-site for the 
progenitor star. It seems likely that the supernova occurred in a stellar wind bubble some 
40kyr ago. The X-ray observations of the pulsar show that the thermal point source has 
a high temperature consistent with the 40kyr age and imply that it was the source of the 
observed SNR. Our deep Chandra pointing reveals extended emission around the pulsar 
with a very compact symmetric structure observed in the 4-8 kev range. The overall PWN 
flux is broadly consistent with the emission expected from this relatively low E spin-down 
pulsar, but the physical origin of the hard emission in these innermost regions is not clear. 
Statistical arguments suggest that the symmetric emission is an equatorial structure, viewed 
edge-on, so that the pulsar spin and motion are roughly aligned as for many other young 
pulsars. However, further observations are needed to reach definitive conclusions. 



We thank Andrew Lyne for providing the current ephemeris of PSR J0538+2817. The 
initial VLBA observations of PSR J0538+2817 were made as part of the VLBA pulsar 
astrometry project (http://www.astro.cornell.edu/~shami/psrvlb/). The Ha image of 
S147 is based on data obtained as part of the INT Photometric Ha Survey of the Northern 
Galactic Plane: prepared by Albert Zijlstra, University of Manchester and Jonathan Irwin, 
loA Cambridge. The National Radio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agreement by Associated Universities, Inc. 
This work was supported by NASA grant NAGS- 13344 and by CXO grant G05-6058 issued by 
the Chandra X-ray Observatory Center, which is operated by the Smithsonian Astrophysical 
Observatory for and on behalf of the National Aeronautics Space Administration under 
contract NAS8-03060. 
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